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PREFACE 


This  report  describes  a  VHF(MB)  link  performance  model  which  was 
constructed  in  three  phases.  The  most  recent  work.  Phase  III,  is  described  in 
Section  I  of  this  report.  This  description  assumes  that  the  description  of 
Phases  I  and  II  have  been  read.  Phase  II  was  described  in  a  draft  report 
which  contained  the  description  of  Phase  I  in  Appendix  A  of  the  report.  The 
draft  report  is  enclosed  as  Section  II  of  this  paper.  Since  the  draft  report 
of  Phase  III  is  enclosed  as  Section  I  of  this  paper,  a  chronological  history 
of  the  model  can  be  obtained  by  reading  Appendix  A,  Section  II  and  I  in  that 
order. 


The  first  phase  resulted  in  a  description  of  meteor  trail  frequency 
and  time  dependence  based  on  limited  experimental  data.  Meteor  trail  frequency 
scaled  with  radiated  power  as  if  all  trails  were  underdense.  Formulas  for 
probability  of  message  receipt  were  developed  for  messages  sent  in  broadcast 
mode  and  split  in  one  to  eight  packets.  The  second  phase  consisted  of  the 
inclusion  of  overdense  meteor  trail  phenomenology  into  the  model.  This 
inclusion  indicated  a  degradation  in  performance  at  lower  radiated  power  than 
predicted  by  the  earlier  model. 


Phase  III  consisted  of  an  examination  of  the  model  for  agreement  with 
experimental  data  and  physics  models  for  trail  frequency  and  time  dependence. 

It  was  decided  that  the  original  modeling  of  these  quantities  was  too  crude  to 
account  for  the  differences  between  the  experimental  situation  where  data  was 
taken  and  the  configuration  of  STRATSEC.  Therefore  a  sky  integration  model 
was  constructed  which  more  accurately  modeled  the  dependence  of  meteor  trail 
frequency  and  time  dependence  on  the  experimental  situation.  The  model  was 
validated  with  some  experimental  data  but  further  data  comparison  is  needed. 
The  new  model  was  constructed  for  surface  to  surface  links.  Extension  of  the 
model  to  the  case  of  one  or  both  nodes  at  balloon  altitude  is  needed  for 
STRATSEC  requirements.  |  Accession  For 
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1.  INTRODUCTION  AND  SUMMARY 


A  VHF(MB)  link  performance  model  was  previously  constructed1  for 
application  to  STRATSEC.  In  this  paper,  we  have  reexamined  the  model  for 

agreement  with  experimental  data  and  physics  models. 

The  first  part  of  the  reexamination  dealt  with  the  predictions  of 
ground  ranges  between  receivers  wherein  meteor  trail  reflection  was  not 

required  for  link  connectivity.  That  is,  the  ground  wave  was  adequate.  The 
original  model  predicted  ground  wave  adequate  for  ranges  up  to  200  km.  We 
obtained  and  exercised  a  ground  wave  path  loss  computer  program  and  found  that 
200  km  was  a  reasonable  bound  for  systems  with  acceptable  path  losses  of  180- 
190  dB  if  the  link  was  over  the  ocean  or  between  antennas  raised  above  poor 
ground.  However,  if  antennas  were  lowered  from  20m  to  2m  over  poor  ground, 
the  ground  ranges  for  acceptable  loss  was  cut  in  half. 

For  the  cases  of  surface  to  FFB  and  FFB  to  FFB  links,  the  predictions 
of  the  original  model  were  verified  by  the  path  loss  computer  program. 

The  second  part  of  the  reexamination  dealt  with  the  performance 
predictions  of  the  meteor  trail  reflection  phenomenon.  This  was  based  on  two 

parameters,  MQ  and  x ,  which  are  the  rate  of  occurrence  of  reflecting  meteor 

trails  of  a  certain  line  density  and  the  effective  life  time  of  the  trail 
respectively. 

The  values  of  these  parameters  as  a  function  of  ground  range  between 

antennas  were  taken  from  the  proposed  model  described  in  a  final  report  by 

2 

Meteor  Communications  Consultants,  Inc.  for  the  Naval  Ocean  Systems  Center. 
Recently,  a  report"1  submitted  by  John  Ames  of  SRI  International  showed 
different  values  of  these  quantities  than  our  model  used.  His  numbers  were 
based  primarily  on  the  work  of  A.  E.  Spezio^  of  the  Naval  Research  Lab. 

While  studying  the  reasons  for  the  discrepancies,  we  found  that  the 
use  of  single  values  of  the  parameters  as  representing  the  entire  area  where 


1 


meteor  trails  occur  is  too  crude  to  effectively  account  for  the  differences 
between  horizontal  and  vertical  polarization  and  the  parameter  boundary  where 
overdense  trails  are  required. 

We  therefore  constructed  a  computer  program  to  sum  contributions  over 
the  sky  where  useful  meteor  trails  occur.  In  addition  to  providing  better 
accuracy  through  the  actual  space  dependence  of  t ,  Mq  and  the  overdense- 
underdense  boundary,  antenna  patterns  can  be  utilized  for  better  reproduction 
of  differing  experimental  situations  for  the  purpose  of  validation. 

The  next  section  describes  the  computer  program  and  how  it  combines 
with  portions  of  the  previous  model  described  in  Ref.  1  to  become  the  new  link 
performance  model.  Reference  1  is  reproduced  as  Section  II  of  this  document. 
An  understanding  of  the  material  in  Ref.  1  is  required  to  understand  this 
section.  The  third  section  shows  the  new  results  and  compares  them  with  the 
predictions  of  the  previous  model  and  with  some  experimental  data. 

In  summary,  comparison  of  the  previous  model  with  experimental  data 
and  physics  models  shows  that  the  model  predicts  over  optimistic  link 
performance. 

We  determined  that  this  was  due  to  the  assumption  of  a  single 
effective  t  at  path  mid-range  and  the  assumption  that  the  overdense  boundary 
depended  only  on  the  power  factor.  The  model  was  modified  to  integrate  each 
path  of  sky  resulting  in  a  more  realistic  link  performance  prediction  for 
surface  to  surface  links  and  the  capability  of  including  antenna  patterns  and 
polarization. 

However,  some  uncertainty  still  exists  for  predicting  performance 
with  vertical  polarization  and  wide  beam  antennas  because  of  the  lack  of 
experimental  data. 

Further  data  is  required  for  model  calibration.  In  particlar, 
vertical  polarization  with  azimuth-omnidirectional  antennas  over  a  series  of 
ranges  and  powers  would  be  useful. 


The  new  model  is  applicable  to  surface  to  surface  links  only. 
Extension  of  the  model  to  FFB  nodes  will  require  further  analysis  and  further 
experimental  data  using  FFBs  or  other  elevated  platforms. 


2.  THE  NEW  LINK  PERFORMANCE  MODEL 


I 

The  new  model  uses  the  same  formulas  for  probability  and  signal  time 
behavior  as  the  old  model  described  in  Ref.  1.  Reference  1  is  reproduced  as 
I  Section  II  of  this  document.  The  major  differences  are  in  the  fact  that  the 

new  model  performs  the  probability  calculations  for  each  portion  of  the  sky 
where  useful  meteor  trails  occur  and  then  combines  the  probabilities,  and  that 
the  signal  reflecting  properties  and  time  constant  of  the  trail  are  more 
j  stringent  than  before. 

For  a  given  value  of  ground  range  between  antennas,  one  quadrant  of 
area  of  the  sky  at  100  km  altitude  is  divided  into  squares  25  km  on  a  side. 
This  quadrant  is  one  fourth  the  area  bounded  by  1000  km  horizons  of  the  two 
antennas.  The  following  quantities  are  evaluated  at  the  centers  of  the 
squares. 


C0S2$,  where  $  is  half  the  scatter  angle; 

EM,  a  quantity  related  to  the  Eshl eman-Manning5  fraction  of  useful 
trails ; 

S0F,  a  quantity  which  contains  the  geometry-dependent  part  of  the 
minimum  trail  line  density  for  signal  threshold. 

These  quantities  are  evaluated  and  stored  for  ground  ranges  from  200 
to  1800  km  in  200  km  intervals.  They  are  called  by  the  new  link  performance 
model  to  evaluate  probabilities  of  message  receipt  at  one  or  more  of  the 
specified  ground  ranges  for  input  values  of  the  engineering  parameters.  The 
input  parameters  are  unchanged  from  those  required  for  the  old  model  described 
in  Ref.  1. 

The  old  model  stipulated  an  effective  rate  of  occurrence  of  meteor 
trails  with  electron  line  density,  q,  of  qQ  =  8.87  •  10*2  electrons/m.  This 


was  about  one  every  two  minutes  at  mid  ranges  and  less  at  longer  and  shorter 
ranges.  This  came  from  data  from  Ref.  2  which  stated  that  a  system  which 
could  tolerate  180  dB  of  path  loss  had  a  meteor  rate  ranging  from  one  every 
minute  to  one  every  three  minutes.  We  chose  two  every  three  minutes  as  a 
compromise.  The  reference  stated  the  overdense-underdense  boundary  which 
corresponds  to  q  =*  qQ  occurred  for  a  system  with  half  the  power  or  any  other 
change  in  engineering  parameters  which  caused  the  tolerated  path  loss  to  be 
reduced  to  177  dB.  By  underdense  scaling  this  meant  the  test  system  described 
required  a  minimum  line  density  q  =  qQ//2  electrons/m.  Therefore  the  meteor 
rate  at  the  boundary  was  MQ  =  2/(3  min)  x  0.700^1/(2  min). 

Unfortunately,  the  linking  of  meteor  rate  to  the  test  system  which 
had  horizontal  polarization  and  high  gain  antennas  and  applying  that  rate  to  a 
system  which  uses  vertical  polarization  and  antennas  which  are  omnidirectional 
in  azimuth  causes  the  results  to  err  on  the  optimistic  side.  This  fact  was 
not  appreciated  until  the  new  model  was  built  which  demonstrated  the 
difference  in  the  two  situations. 


The  new  model  obtains  the  meteor  rate  as  follows.  Eshleman  and 

5 

Manning-'  state  the  number  of  meteors  per  unit  area  per  second  entering  the 
earth's  atmosphere  with  line  density  q  or  higher  is 

2 

N  =  C/q  meteors/km  /sec 

O 

where  C  ^1.6*10  ,  a  number  we  can  use  but  also  decrease  to  account  for  times 
of  day  or  year  when  fewer  are  expected  in  the  region  of  interest.  Of  these,  a 
fraction,  p,  are  oriented  in  such  a  way  as  to  produce  a  reflection.  The 
formula^  is 


50  km 

P  x 


EM 


5 


where  h 


=  meteor  altitude  =  100  km 


D  =  1/2  ground  range,  R  (RH2D) 

S  -  (Rt  +  Rp)/R 

n  «  (Rt  -  RP)/R 

Rt>Rp  -  distances  from  antennas  to  point  in  sky. 

Since  there  are  four  patches  of  sky  with  the  same  geometry,  the  total 

o 

number  of  meteors/min  of  useful  orientation  for  a  (25  km)  patch  of  the 

2 

quadrant  with  a  value  of  q  =qQ  *  (4pC/qQ)  (25  km)  (60  sec/min).  Therefore  Mq 
=  1.8*10'8C  (EM)/R  meteors/min  is  the  required  rate. 


The  second  quantity  of  concern  is  the  meteor  burst  time  constant. 


2  2 
A^  sec  fl 

16-rr2d 


(sec) 


where  A  =  wave  length  (m) 


d  =  diffusion  coefficient  (itr/sec), 


This  is  evaluated  in  the  new  model  for  each  patch  of  sky.  As  suggested  in 

Ref.  2,  the  old  model  used  a  single  value  of  x  where  $  was  obtained  for  a 

point  at  100  km  altitude  midway  between  the  antennas.  Although  this  point  has 

the  largest  value  of  sec  $  and  hence  x  of  the  sky,  it  was  believed  at  the 

time  that  most  meteor  trails  of  importance  occurred  near  that  point.  However, 

the  new  model  shows  many  contributing  meteor  trails  over  a  large  region  of  sky 

2 

which  brings  down  the  effective  value  of  sec  $  and  hence  x. 


The  new  model  also  shows  the  difficulty  of  interpreting  experimental 


data  on  t  for  the  purpose  of  fixing  the  value  of  d.  We  have  used  d  =  8m  /s  in 

the  past.  The  appropriate  value  depends  on  the  altitude  of  the  trail.  The 

2 

diffusion  coefficient  varies  from  0.5  at  80  km  to  100  m  /s  at  100  km  altitude. 

An  experimental  determination  of  x  by  Keary  and  Wirth6  show  a  two  to  one  daily 

variation  of  x  which  they  attribute  to  higher  velocity  meteors  in  the  early 

morning  causing  the  trails  to  occur  higher  in  the  atmosphere  where  d  is 

greater.  Their  experiment  used  horizontal  polarization  and  Yagi  antennas 

which  lead  to  slightly  higher  values  of  effective  x  than  would  obtain  for 

vertical  polarization  and  omnidirectional  antennas.  We  illustrate  this 

difference  by  using  the  code  to  calculate  an  effective  x  by  averaging  over 

the  sky  weighted  by  EM/S0F  which  represents  the  relative  number  of  useful 

underdense  trails  for  a  specific  system.  We  show  in  Figure  1,  the  results  of 

2 

the  calculation  for  a  frequency  of  47  MHz  and  d  =  8m  /sec.  The  circles  are 
the  old  midpoint  calculation  whose  values  are  too  high.  The  B's  are  the 
result  of  simulating  an  experiment  by  Boeing  described  in  Section  3  which  used 
horizontal  polarizations  and  Yagi  antennas.  The  x's  result  from  vertical  stub 

3 

antennas  which  cover  a  larger  region  of  sky.  The  dots  are  from  Ames'  report  . 

The  Keary  and  Wirth  data  bar  at  700  km  range  shows  agreement  at  the  low  end 

2 

but  suggests  that  we  could  reduce  d  to  4m  /sec  for  early  morning.  One  other 

data  point  from  the  COMET  system  (described  in  Section  3)  at  1000  km  falls 

2 

directly  on  the  B  so  we  will  stay  with  d  ■  8m  /sec. 

The  final  quantity  used  for  the  probability  calculations  is  SQ.  This 
is  the  quantity  defined  in  the  old  model  as  a  required  signal  as  determined  by 
system  parameters  normalized  to  the  overdense-underdense  boundary  at  qQ.  For 
$0  <1,  underdense  trails  sufficed.  In  the  old  model,  SQ  was  matched  to  the 
data  of  Ref.  2.  In  the  new  model,  SQ  is  taken  directly  from  the  radar 
transmission  equation  which  includes  the  meteor  trail  scattering  cross  section. 

PtGtGr  16Tr2  RtRr(Rt  +  V  (1  ‘  c°s2Bsin2$) 

Tr  xVr2  sin2a  gt(0t.<J>t)gr(9r»<frr) 


The  quantity  on  the  left  is  the  basic  path  loss  (when  in  dB)  that  the 
system  can  tolerate  to  achieve  the  required  threshold  power,  T  .  The  radiated 
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o 
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power  is  and  Gt  and  G  are  the  maximum  gains  of  the  antennas  over  an  isotrope. 

The  quantity  on  the  right  is  that  loss  in  scattering  from  the 
properly  oriented  meteor  trail  of  line  density  q  where  gt  and  gp  are  the 
normalized  antenna  gain  patterns  of  transmitter  and  receiver.  These  are  to  be 
evaluated  for  the  elevation  and  azimuth  angles,  9  and  4>,  to  the  point  in  the 
sky.  The  angle  a  is  between  the  electric  vector  at  the  meteor  trail  and  Rp, 
the  direction  vector  from  meteor  trail  to  receiver  and  B  is  the  angle  between 
the  trail  and  the  plane  of  Rt  and  Rp. 


We  choose  the  value  of  SQ  as  the  normalized  line  density  required  for 


the  experimental  and  geometrical  parameters  involved. 

V  -  Tr"  T  SOF(S) 

"  n  n  r 


IS. 

0  % 


w 


where 


S0F(3) 


.  £ 


+  Rr)(l  -  cos26sin2$>) 


— * 

sin  a  gt(et,<i>t)gr(er,<t>r) 


by  substitution  from  the  transmission  equation. 

Since  Eshleman  and  Manning5  integrated  over  B  to  obtain  their 
fraction  of  properly  oriented  trails,  we  must  integrate  over  B  to  obtain  the 
proper  average  for  the  probability  calculations  without  knowing  how  to  weight 
the  values  of  B  in  the  averaging  process.  The  proper  average  is  of  the 
quantity  1/S0F(B)  or 

1 

VI  -  cos2B  sin2$ 

which,  for  trails  uniformly  distributed  in  B,  has  the  average 
^  K  (sin  $>) 


where  K  is  the  complete  elliptic  integral  of  the  first  kind.  We  choose 

instead  the  simpler  quantity 

2 

1  +  cos  $ 

which  is  close  to  the  average  but  smaller  for  large  $  making  the  model 
conservative.  The  true  weighting  over  6  is  unknown  without  redoing  the  work 
of  Eshleman  and  Manning.  We  therefore  use 


SOF  = 


_  (i  +  cos 
~  2  sin  a 


WRt  +  Rr> 


g(et  ,d>t  )g(er,q>r  j 


For  vertical  stick  antennas. 


g(9,<J>)  =  cos^e 


and  the  maximum  gain  over  an  isotrope,  G  *  3/2. 


Finally,  the  quantities  SQ,  MQ  and  t  as  specified  are  used  as 
described  in  Ref.  1  to  calculate  the  probability  of  not  receiving  a  signal  of 
required  value  SQ  and  time  length  At  (no  useful  meteor  bursts)  assuming  no 
useful  meteor  bursts  occurred  prior  to  a  time  T  before  the  signal  was  sent. 
This  is  done  for  each  part  of  the  sky  and  their  product  is  then  the 
appropriately  combined  probability  for  the  entire  sky.  This  is  used  to  obtain 
the  probability  of  message  receipt  as  described  in  Ref.  1. 


3.  PREDICTIONS  AND  DATA  COMPARISON 


The  new  model  has  the  capability  of  being  adjusted  to  reproduce 
different  experimental  situations.  We  illustrate  this  by  comparing  it  with 
some  performance  measurements  taken  by  Boeing  described  In  NOSC  TR1387.  Data 
was  taken  for  two  paths  of  800  and  1600  km  ground  range.  A  probe  signal  was 
continuously  sent  from  a  master  to  a  remote  station.  If  a  meteor  trail 
occurred  of  sufficient  strength  for  the  remote  station  to  detect  the  probe 
signal,  the  remote  station  sent  a  small  block  of  data  to  the  master  station 
over  the  same  trail.  The  performance  data  was  in  the  form  of  probability  of 
success  vs.  waiting  time.  This  data  was  averaged  over  a  two  week  period  for 
six  4  hour  time  blocks  to  show  the  daily  variation  of  meteor  rate. 

The  allowable  path  loss  for  the  probe  signal  was  191  dB.  The 
allowable  path  loss  for  the  response  signal  was  199.7  dB.  The  probe  signal 
required  a  stronger  meteor  trail  by  8.7  dB  primarily  because  of  the  noisy 
remote  receiver. 

This  means  the  occurrence  of  a  meteor  trail  strong  enough  for  the 
remote  receiver  to  detect  the  probe  signal  will  allow  the  remote  station  to 
send  the  short  data  block  (which  required  At  *  0.156  seconds)  on  the  same 
trail  whose  scattering  cross  section  would  have  diminished  by  a  factor 
e'2 


For  8.7  dB,  t  ^0.156  seconds  is  sufficient  for  the  correct  message 
receipt  to  be  determined  solely  by  the  probe  path  loss  of  191  dB.  From  Figure 
1,  the  Boeing  values  of  t  are  greater  than  0.2  sec. 

In  the  expression  for  S0F,  on  page  10  we  alter  the  calculation  of 
sin  a  to  account  for  horizontal  polarization  and  use  the  experimental  values 
of  Gt  *  20  (13  dB)  and  Gp  =  6.3  (8  dB)  in  an  approximate  expression  for  the 
normalized  gain  pattern 

e-G(92  +  <J>2)/4 
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The  first  data  check  on  the  new  model  is  the  ratio  of  average  waiting 
time  for  the  two  ranges  800  and  1600  km.  The  data  show  the  waiting  time  ratio 
to  be  about  2  (the  longer  range  has  the  larger  waiting  time).  This  is  in 
agreement  with  the  model  for  a  191  dB  system  but  not  lower  path  loss  systems. 
The  old  model  also  showed  a  factor  of  2  but  had  no  dependence  of  the  factor  on 
allowable  path  loss.  Since  the  output  can  be  put  in  the  form  of  useful  meteor 
rate  M  through  the  formula  for  probability  of  success  vs.  waiting  time  t, 

pal-  e-Mt  (1 ) 


we  match  their  curves 

values  of  M: 

Of  P 

vs.  t 

to  the 

formula  and  obtain 

the 

Time  of  Day 

0-4 

4-8 

8-12 

12-16  16-20 

20-24 

M 

1.6 

2.3 

1.7 

1.0  0.75 

0.6 

C(xl08) 

0.9 

1.3 

0.94 

0.55  0.4 

0.33 

The  third  row  contains  the  corresponding  values  of  C  by  which  the 
model  would  predict  the  values  of  M.  We  recall  that  C  was  associated  with  the 
overall  rate  of  meteors  entering  the  earth's  atmosphere  which  Eshleman  and 

C  O 

Manning  estimated  as  1.6*10  .  The  experimentally  determined  values  of  C  are 
slightly  lower  but  very  close.  The  old  model  would  have  predicted 

M  =  Mq  10  *05  (U  dB)  =  2/3  (3.55)  -  2.4 
which  is  too  high. 

The  experiment  was  In  June  and  showed  a  daily  variation  of  meteor 
rate  of  4  to  1.  July  and  August  are  supposed  to  have  the  highest  rates  and 
February,  the  lowest.  The  ratio  of  summer/winter  rates  is  not  precise, 
ranging  from  1.4  in  Ref.  2  to  5  in  Ref.  7. 


In  order  to  demonstrate  the  difference  between  the  old  and  new  model, 
we  reproduce  a  curve  of  probability  of  correct  message  receipt  vs.  power 
factor  P^G^Gp  from  Figure  14  of  Ref.  1  in  Figure  2.  Using  the  identical 
parameters  in  the  new  model  for  short  stub  vertical  antennas  and  a  value  of 

Q 

C  =  1.4  10  as  the  incident  meteor  rate  at  the  best  time  of  day  and  year  we 
obtain  the  curves  shown  in  Figure  3.  The  new  model  shows  much  higher  power 
factors  required  throughout.  Some  of  the  difference  is  due  to  the  more 
conservative  values  of  MQ  and  x,  but  a  major  difference  is  the  new  method  of 
treating  the  underdense/overdense  boundary.  The  old  model  put  the  boundary  at 
a  500  watt  power  factor.  This  is  the  "waterfall"  region  of  Figure  2.  The  new 
model  shows  overdense  bursts  to  be  required  everywhere  for  a  1400  km  range 
until  the  power  factor  becomes  greater  than  2500  watts.  We  illustrate  this 
point  by  contour  plots  of  SQ  vs.  sky  position  for  three  ground  ranges  of  600, 
1000  and  1400  km  in  Figures  4,  5  and  6  for  a  180  dB  system  at  47  MHz.  In  the 
examples  of  Figures  2  and  3  this  corresponds  to  a  power  factor  of  1000  watts. 
The  coordinates  are  distances  from  the  midpoint  in  km.  The  abscissa  is  along 
the  line  between  antennas.  In  Figures  4  and  5,  the  region  above  one  antenna 
is  shown  to  require  high  values  of  SQ  due  to  polarization  losses.  This 
doesn't  happen  for  horizontal  polarization.  The  overdense  boundary  is  at  SQ  * 
1.  Since  SQ  scales  as  the  square  root  of  power,  at  4000  watts  the  boundary 
would  be  at  the  contour  marked  SQ  =  2.  According  to  Figure  3,  this  power 
increase  is  sufficient  to  increase  all  probabilities  to  greater  than  0.96. 

We  have  used  the  best  meteor  rate  for  a  comparison  with  the  old  model 
and  found  the  old  model  too  optimistic.  We  will  choose  a  worst  rate  as  one 
fourth  the  best  to  typify  a  bad  time  of  day  in  a  good  month  or  vice  versa. 
The  results  are  shown  in  Figure  7.  Now  much  higher  powers  are  required  for 
the  same  broadcast  time. 

Since  all  results  are  functions  of  the  product  of  overall  meteor 
rate,  C,  and  broadcast  time,  t.  We  can  reproduce  Figure  3  for  the  worst  rate 
by  increasing  broadcast  time  to  40  min. 

Finally,  we  illustrate  the  best  case  code  predictions  for  a  single 
packet  of  512  bits  at  a  data  rate  of  4000  bps  at  47  MHz  for  the  vertical  stub 


400  KM.TBR 


Probability  of  Correct  Message  Receipt  Vs.  Power  Factor  P.G.G  ;  New  model.  Worst 


case  in  the  form  My  vs.  range  in  Figure  8.  The  various  curves  show  the 
performance  gain  obtained  by  increasing  the  allowable  path  loss  Lg.  For  the 
case  considered,  an  Lg  of  180  dB  translates  to  a  power  factor  of  PtGtGp  a  1000 
watts.  Therefore  190  dB  corresponds  to  10,000  watts.  These  curves  can  be 
used  to  obtain  the  probability  of  message  receipt  for  a  broadcast  time  t  by 
expression  (1)  on  page  12.  For  example  at  1400  km,  the  180  dB  curve  shows  Mv 
=  0.112/min.  Therefore  a  broadcast  time  of  10  min  produces  a  probability  p  = 

i  ip 

l-e  *  =  0.67  as  seen  on  Figure  3.  To  obtain  the  worst  case,  we  divide  M 

0  28  V 
by  4  to  obtain  p  =  l-e  =  0.244  as  shown  on  Figure  7. 

The  curves  of  Figure  8  indicate  that  higher  power  factors  than 
previously  thought  are  needed  to  overcome  the  overdense  effects  at  the  longer 
ranges. 

Although  the  results  presented  here  are  pessimistic  compared  to  the 
old  model,  they  should  get  better  at  long  ranges  by  the  use  of  antennas  with 
better  gain  patterns  than  the  stubs  used  here  for  illustration.  Also,  the 
effect  of  one  or  more  elevated  nodes  has  not  been  analyzed.  The  old  model 
assumed  the  effect  to  be  the  same  maximum  meteor  rate  at  midrange  with  a 
longer  extent  until  long  range  cut  off.  The  mid  range  could  actually  turn  out 
to  have  a  higher  meteor  rate  because  of  the  longer  horizon  and  the  fact  that 
portions  of  the  sky  far  from  the  midpoint  are  useful  for  vertical 
polarization. 
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1 .  INTRODUCTION 


A  VHF(MB)  link  performance  model  was  previously  constructed^  for  appli¬ 
cation  to  STRATSEC  based  on  the  physics  of  underdense  reflection  of  radio  signals 
by  meteor  trails.  This  report  describes  the  improvement  of  that  model  by  the 
inclusion  of  overdense  reflection  characteristics.  The  model  improvement  removes 
erroneous  high  predictions  of  performance  by  non  existent  underdense  trails  which 
the  model  required  for  low  power  levels. 

All  underdense  models  assume  the  signal  reflected  by  the  trail  to  be 
proportional  to  q,  the  electron  line  density. 

S  =  reqe‘t/T  (1) 

-15 

where  r  =  2.8  •  10  m  =  classical  electron  radius  and  t  is  the  trail  time 

e  2 
constant.  The  overdense  variation  of  the  signal  is 


S  = 


(4ir)2x 


Jin 


'4reqx 


which  has  a  maximum  value  of 


and  drops  to  zero  at  t  =  4rgqx  when  the  trail  is  completely  underdense. 

Any  meteor  trail  can  be  considered  to  have  an  overdense  core  and  an 
underdense  outer  shell  at  early  time  with  the  overdense  core  shrinking  to  zero 
at  t  =  4rgqx.  For  small  values  of  4rgq,  the  overdense  signal  is  short  lived  and 
the  dominant  behavior  is  given  by  Eq.  (1).  The  previous  model  assumed  this 
variation  for  all  values  of  q  and  t. 

For  large  values  of  4rgq,  however,  the  overdense  core  is  the  actual 

reflection  mechanism  until  t  =  4r  qx  >  x  after  which  some  residual  underdense 

e 

signal  is  possible.  In  this  case  the  overdense  signal  maximum  value,  $m,  is  less 
than  rgq  and  the  real  signal  value  is  smaller  than  the  underdense-only  theory 
prediction. 
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This  impacts  the  link  performance  model  in  the  following  way.  The 
output  of  the  model  is  the  probability  of  receiving  a  message  by  meteor  burst 
reflections  which  depends  on  the  quantity  M,  the  average  rate  of  useful  meteor 
trails.  Since  the  frequency  of  trails  is  inversely  proportional  to  the  square 
of  their  line  density,  the  higher  a  line  density  required  to  reflect  a  signal, 
the  lower  the  value  of  M. 

The  original  model  obtained  the  value  of  M  by  scaling  from  data  obtained 

3 

from  underdense  trails.  From  that  data,  Mq0,  the  average  rate  of  useful  meteor 
trails  was  determined.  Reference  3  ''••tlined  a  method  of  scaling  from  their  test 
system  to  any  other  system  operating  at  the  same  optimum  ground  ranges.  For  a 
required  receiver  threshold  Tr(watts)  which  sets  the  required  received  signal 
amplitude  /T^  and  given  values  of  transmitter  gain,  Gt,  above  an  isotrope,  trans¬ 
mitter  radiated  power,  Pt;  receiver  gain,  Gr»  above  an  isotrope  and  wavelength,  X, 
the  required  value  of  signal  reflected  by  the  trail  (*  /cross-section)  is  propor¬ 
tional  to 


3 

For  the  test  system  , 


Therefore,  the  required  value  of  signal  for  any  other  system  scales  as 


For  underdense  bursts,  $0Q  *  q  and  since  the  average  rate  of  meteors 
with  line  density  q  or  higher  is  proportional  to  1/q,  the  scaling  for  underdense 
bursts  is 
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For  a  system  with  a  large  value  of  SQ0  where  the  corresponding  value  of  q  is  so 

large  that  only  the  overdense  reflection  is  operative,  the  signal  is  «  and  not 

4 

q.  In  this  case  the  required  q  scales  as  SQQ  and  the  rate  of  useful  meteors  is 
M  «  which  is  lower  than  the  underdense  scaling.  With  lower  values  of  M,  the 
probability  of  message  receipt  will  be  lower  than  predicted  by  the  underdense 
scaling. 

The  next  section  describes  the  analysis  of  overdense  theory  and  the 
improved  link  performance  model.  The  third  section  shows  the  new  results  and 
compares  them  with  the  predictions  of  the  old  model.  The  description  of  the 
old  model  in  Ref.  1  is  included  as  Appendix  A  and  Appendix  B  describes  the  imple¬ 
mentation  of  the  new  model  as  a  computer  program. 
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2.  ANALYSIS 


In  this  section,  the  transition  point  for  underdense-overdense  reflections 
is  obtained  and  the  method  of  describing  overdense  signals  and  the  impact  on  prob¬ 
ability  calculations  is  described. 

Reference  3  described  the  transition  point  to  be  at 
M/Mqo  *  0.7  . 

This  value  agrees  with  our  conjecture  in  our  description  of  overdense  effects  in 

reference  1  that  the  transition  would  occur  at  S_  ^  /T. 

oo 

In  our  previous  engineering  model1,  the  value  of  MQQ  was  constant  over 
most  of  the  range  with  a  linear  fall  off  at  large  range  due  to  shadowing  of  the 
useful  area  of  sky  by  the  curvature  of  the  earth.  Since  this  shadowing  has  nothing 
to  do  with  required  signal  strength,  the  overdense  scaling  should  apply  to  the 
long-range  diminished  MQQ  also. 

In  addition,  there  was  a  slight  fall  off  at  short  ranges  for  surface  to 
surface  links  only.  This  fall  off  was  attributed  to  geometric  and  polarization 
effects.  This  means  the  required  signal  strength  is  greater  at  the  short  ranges 
and  the  scaling  should  be  different.  We  have  insufficient  information  as  to  all 
the  reasons  for  the  short  range  fall -off  to  know  how  to  change  the  scaling.  Until 
we  obtain  that  information  we  will  apply  the  overdense  scaling  to  the  short-range 
diminished  M  also. 

We  define  a  dimensionless  q  which  is  unity  at  the  transition  point  SQ0  *  /T. 

We  scale  to  the  transition  point  by  defining  *  M L  //F,  «  S so  that  the 

3  o  oo  o  oo 

average  rate  of  meteors  with  line  density  q  or  higher  is  M  =  MQ/q  and  the  signal 
reflected  by  a  trail  of  line  density  q  has  a  maximum 

S(q)  =  q  q<l 

S(q)  *  q*  q>l  . 

The  required  signal  as  determined  by  the  system  parameters  is  SQ.  For 

S  <<1 ,  underdense  meteor  trails  suffice.  For  S  >1 ,  overdense  meteor  trails  are 
o  o— 

required.  If  SQ  is  near  and  below  unity,  overdense  trails  may  still  be  required. 


Since  we  wish  an  entire  message  or  packet  of  length  At  to  be  reflected  by  a  single 
trail,  the  required  reflecting  signal  must  be  greater  than  SQeAt//T  since  an  under- 
dense  trail  has  reflecting  signal  time  dependence 

S  =  qe-t/T  . 

The  overdense  signal  has  a  complicated  time  dependence  which  is  obtained 
from  an  approximate  theoretical  treatment.  The  important  characteristics  are  its 
maximum  and  the  time  at  which  it  goes  underdense  which  is  proportional  to  q. 

Examples  of  signals  from  overdense  bursts  are  shown  in  Figure  1,  an 

extract  from  a  chart  recording  obtained  at  a  receivino  site  over  a  60  mile  path 

4 

using  a  10  watt  CW  transmitter.  The  large  seven-second  signal  and  the  final  two 
are  overdense.  The  many  small  and  one  larger  signal  following  the  seven  second 
signal  are  underdense.  The  fading  in  the  large  signal  makes  it  less  useful  than 
one  would  predict  from  its  average  envelope  value.  It  appears  that  short  signals. 
At  <  1  sec  have  a  high  probability  of  transmission  if  the  required  strength  is  half 
the  envelope  value.  The  small  overdense  signals  do  not  exhibit  this  fading  and  in 
fact  appear  closer  to  the  exponential  time  behavior  of  underdense  signals. 

In  order  to  be  conservative  and  for  ease  of  computation,  we  will  treat 
the  overdense  signals  as  if  they  also  had  an  exponential  time  decay  but  with  a 
longer  time  constant  associated  with  the  overdense  period. 

The  transition  will  be  smooth  if  we  choose  the  time  behavior  of  the  over- 
dense  signal  as 


t  _<  -rq 


h  q-1  -t/x  * 

e  '  t  > 


which  drops  to  1/e  of  its  initial  value  at  the  end  of  the  overdense  period  and  then 
diminishes  further  as  an  underdense  trail. 

We  now  outline  the  procedure  for  the  determination  of  the  probabilities 
of  message  receipt.  As  shown  in  the  original  model  description^,  all  possible 
probabilities  of  interest  are  obtainable  from  the  following  probability: 


5 


Figure  1.  Extract  from  Chart  Recording  at  Receiving  Station  Showing  Overdense  and  Underdense  Type  Signals 


F(SQ,At,T)  is  the  probability  of  not  receiving  a  signal  of  required 
value  SQ  and  time  length  At  (no  useful  meteor  bursts)  assuming  no  useful  meteor 
bursts  occurred  prior  to  a  time  T  before  the  signal  was  sent. 

Because  of  the  change  in  the  assumed  time-dependence  of  the  signal 
strength  depending  on  whether  the  meteor  trail  is  overdense  or  underdense,  several 
cases  need  to  be  considered  in  order  to  calculate  the  relevant  probabilities. 

•t/x 

These  cases  depend  on  whether  the  amplitude  of  a  reflected  signal  decays  as  e 
(an  underdense  trail),  e-t^qx  (an  overdense  trail),  or  changes  character  during 
the  relevant  time.  Probabilities  for  the  different  cases  are  described  below.  The 
situation  is  diagramed  in  Figure  2  for  the  underdense  case  where  SQeAt//T  «  1. 


Figure  2.  T  <  tQ,  Underdense  Only. 


Because  of  the  time  dependence  of  the  reflected  signal,  the  minimum  q 
required  for  sufficiently  reflecting  the  signal  depends  on  the  time  frame  prior 
to  the  time  the  signal  was  sent.  The  time  t  corresponds  to  the  prior  time  which 
requires  a  minimum  value  of  q=l ,  the  transition  value, 

t  ■  -  t  Jln(SQ)  -  At,  a  positive  number  . 

The  fact  that  T  <  t  means  we  are  not  concerned  with  overdense  bursts  since  we  are 

o 

given  that  no  useful  bursts  occurred  prior  to  T  before  the  signal.  In  this  case 
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the  formula  for  the  probability  is  that  derived  in  Ref.  1.  We  will  rederive  it 
here  to  illustrate  the  procedure. 

For  a  meteor  with  line  density  in  an  interval  dq  about  q,  the  average 
2 

rate  is  (MQdq/q  )  and  the  prior  time  interval  within  which  this  meteor  is  useful 
is 

t(q)  =  t  in  (q/qQ) 


where 


V 


At/x 


is  the  minimum  useful  value  of  q.  For  the  given  time  and  the  given  average  rate, 
the  probability  that  a  meteor  of  that  type  does  not  occur  in  that  time  is 


P(q) 


t(q) 


We  assume  the  occurrence  of  a  meteor  of  one  value  of  q  to  be  independent 
of  the  occurrence  of  a  meteor  with  another  value  of  q  and  multiply  the  probabilities 
of  nonoccurrence  of  the  various  types. 

For  values  of  q  ^  qM  =  qQeT/,T  such  that  t(q)  >_  T,  the  probability  of  non¬ 
occurrence  is 


P(q) 


e 


T 


since  none  occurred  prior  to  T  before  the  signal. 

The  total  probability  of  nonoccurrence 


pu(S0,At,T)  =  e  \  q 


/*qM  M  dq  f°° 

-  f  x*n  (q/q  )  -J qM 

=  p  •'n  n  P 


M  dq 


=  e 


Hi  t/ 

-  f-  0-e'T/T 


) 


as  derived  previously  in  Ref.  1. 
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V  ^ \  ■  .  • 


The  fact  that  we  have  integrated  beyond  q=l  means  that  we  have  included 
overdense  bursts  in  the  probability  calculation.  They  have  been  correctly  included 
since  all  we  need  to  know  is  that  an  overdense  burst  occurring  within  time  T  of 
the  signal  will  be  useful  and  that  is  automatically  true  for  T  <  tQ. 

The  underdense  reflection  situation  for  T  >  tQ  is  diagramed  in  Figure  3. 


Figure  3.  T  >  tQ,  Underdense,  Plus  Overdense. 


For  this  case,  the  probability  of  no  useful  meteor  is  the  product  of 
the  probability  of  no  useful  meteor  in  time  tQ  (which  is  an  underdense  situation) 
and  the  probability  of  no  useful  meteor  of  q  >  1  in  the  remaining  time  (T-tQ).  A 
glance  at  the  figure  shows  the  latter  probability  to  be  identical  to  the  prob¬ 
ability  that  would  be  calculated  in  an  overdense  case  for  prior  time  (T-tQ)  and  a 
signal  of  length  At+tQ.  Denoting  F  for  the  calculated  probabilities  where  qQ  >_  1 
(overdense  case),  we  have  F  =  F^(S0,At,tQ)  F0(S0>At+t0,  T-t  ). 

In  the  general  case  of  overdense  trails,  we  will  have  to  include  the 
probability  of  occurrence  of  meteors  with  q  >_  q^  where  q-j  ^1  is  the  smallest 
relevant  meteor  trail  strength.  (In  the  above  case  we  know  q-j  =  1.)  For  any  q  1, 
the  signal  form 
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or  until  t(q)=T  which  occurs  for 


(T+At) 


qT 


%  e  qTx  =  S, 


For  the  case  qT<q2  or  <q2*  we  have 


f  -q_  M  dq  /  /A\  \  /*“  H  dq 

-  I  qTAn  s-  -At)  +  i  -VT 

P^(S0,At,T)  =  e  [\  q2  V  V V  )  \  q 


Mr  /  0  o  \  /qT\  MnAt 

=  g-  —  Un2(qT)-«n2(qi)j  jm(SQ)  in  ^(T+AtJ/q., 


For  A  =  and  B  = 

qyt  q-|T 


and  the  definitions  of  q-j  and  qT  this  can  be  put  in  the  form 


_  -M  t (A-B ) ( 1 +2[A+B] ) 


P  =  e 
*0  e 


For  the  case  ^2 

qT>q2  and  we  have  the  same  integration  until  q  =  q2.  For  q>q2  we  use 


q%  eq-le-(t+At)/r  .  s 


or  t 


/  / fM \\  qM-i  -(T+4t)/T 

=  ilq-1  +  An  If-)]-  At  until  qM  which  occur  forqM4e  e  =  S 


so  that 


-M  T 
r-  e  0 


£n2(q2)-in2(q1) 
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£n  (SQ)  £n 


(qi)  **  fc 


-vrM  4 

*o  0  •'^2  "2 


q-l  +  *n  (Pi  ^ 

Vo/  T 


The  integral  is  £n 


4lt!-M .  4o 

Since  qM  *  q2e  T  , 


*"‘e<U  *"<eq2>  ,  /T+4t 


4<! 


m 


and  by  the  definition  of  q2  and  q^ 


-Mr 

0 


3-B-2B2  *  4qm(c-l )  ♦  } 


V>2  %j 


where 


r  _  T+At 

"  V  ’ 


Finally  for  At/x>q2,  we  have  the  situation  of  Figure  5  where 

t(q)  *  x|q-l  +  £n  -  At /x^ 
and  the  integration  is  between 

q3  =  q2e4^At/T"q3^  and  qM 


__  -Hr 
and  =  e  0 


£n 


■M0t 


£n(5n 


1  )  1 


q3/  4  \qm  q: 
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)  in(eqM) 


tn(eqw)\  , 
%,  / 


which  reduces  to 


3.  RESULTS 


The  implementation  of  the  overdense  formulas  is  described  in  Appendix  B. 
We  now  illustrate  the  results  of  the  new  model  for  cases  of  interest  previously 
described  in  Ref.  1. 

Figure  6  is  an  illustration  of  a  FFB-FFB  link  PMR  vs.  broadcast  time 
for  a  value  PtGtGr  =  10°  watts.  Two  of  the  curves  are  the  N=1  and  N*8  packet 
cases  of  the  old  model.  The  third  and  lowest  curve  is  the  N=8  packet  case  treated 
by  the  new  model.  The  N=1  case  output  is  identical  for  the  two  models.  We  see 
the  difference  as  negligible  for  this  case  since  the  meteor  trails  of  importance 
are  primarily  underdense. 

In  that  reference  we  suggested  a  scaling  of  power  requirements  with 

2 

broadcast  time  as  P^G^  t  =  constant  which  comes  from  the  theory  of  underdense 
trails.  We  therefore  suggested  that  results  for  1000  watts  are  identical  to 
results  for  250  watts  if  the  time  is  doubled.  We  now  show  the  underdense  scaling 
to  overstate  performance  when  carried  to  powers  as  low  as  250  watts.  This  results 
because  the  required  trails  are  overdense  at  that  power.  Figure  7  shows  the  new 
results.  The  top  two  curves  show  the  erroneously  scaled  results  for  one  and  eight 
packets.  These  are  the  same  curves  as  in  Figure  6  but  the  time  scale  (abscissa) 
has  been  doubled.  For  N=1 ,  P=0.9  at  t=8  min  according  to  the  old  scaling. 

The  lower  two  curves  are  the  actual  results  from  the  new  model.  The 
time  for  P=0.9  is  much  longer.  For  the  N=1  case  we  can  explain  the  difference. 

The  probabilities  depend  on  the  product  Mt  where  M,  the  useful  burst  rate,  =  MQ/q 
where  q  is  the  required  strength.  For  1000  watts,  we  had  q  so  M  =  MQ/2". 

The  product  P  =  0.9  was  Mq  /?  x  4  min.  For  250  watts  the  required  signal  level 
was  twice  as  high  which  would  require  q  *  if  underdense  scaling  were  valid 
leading  to  a  product  (MQ//?)  t  *  Mfl  /?  x  4  or  t=8  min  for  250  watts.  The  new  pro¬ 
gram  recognized  that  overdense  bursts  are  required  for  signal  levels  >1  so  the 
required  strength  is  q4  =  /2  or  q=4  leading  to  M  =  MQ/4.  The  required  product 
(M(/4)t  =  Mq/?  x  4  or  t  =  16  /S"  -  22.6  min  as  seen  in  Figure  7.  This  scaling 
argument  may  hold  for  N*1  since  the  time-dependent  form  of  the  trail  signal  is 
not  important  if  At«x  which  is  the  case  here.  For  N=8,  8At  %r  and  the  different 
time  dependences  of  over-  and  underdense  trails  complicate  the  problem  and  the 


500  KM 


Figure  7.  Probability  of  Correct  Message  Receipt  Vs.  Broadcast  Time  FFB-FFB,  PtGtG, 


scaling  argument  becomes  inaccurate.  This  new  scaling  argument  is  also  inaccurate 
for  any  N  for  the  break  point  case  of  500  w.  For  this  case  the  underdense  scaling 
argument  would  say  the  required  strength  is  q=l  leading  to  t  *  4  /7  *  5.66  min. 

The  N*1  and  8  cases  are  shown  in  Figure  8  for  500  watts.  The  actual  time  to  achieve 
P  =  0.9  with  N=1  is  approximately  7  min  when  overdense  trails  are  considered. 

The  reason  for  this  discrepancy  is  that  At/x  is  finite.  For  finite 
At/x,  the  minimum  strength  required  is 

q  =  SoeAt/T  q  <  1 
4At 

q  =  SQ4e  qx  q  >  1 
For  At/x  s  0.1  we  have 

q(1000)  =  -^e0'1 

n  /z 

4At  4At 

q(250)  *  (/T)4  eqT  Jfc  4  e"5*"  %  4  e0,1  £4 

4At  0^ 

so  the  scaling  that  ignores  At/x  is  adequate.  But  q(500)  =  eqT  *  e  q  is  not 
accurately  approximated  by  q=l .  The  answer  is  q  *  1.346. 

We  therefore  abandon  scaling  arguments  and  use  output  curves  of  prob¬ 
ability  of  message  receipt  vs.  power  for  a  five  minute  and  ten  minute  broadcast 
time.  Figures  9-14  show  the  results  for  the  two  times  and  the  three  cases  of 
FFB-FFB  at  2500  km,  FFB-GRD,  and  GRD-GRD  at  1400  km.  These  three  cases  were  repre¬ 
sented  in  the  power  requirements  section  of  Ref.  1.  We  replace  the  values  of 
radiated  power  given  in  Table  5-2  of  that  reference  with  the  values  obtained  from 
the  N=7  curve  in  Figures  9-14.  In  these  figures,  the  abscissa  is  the  power  factor 

P..GJL  for  the  FFB-FFB  and  GRD-GRO  cases  but  is  PJa+  for  the  FFB-GRD  case  where  an 
t  t  r  t  t 

average  over  range  was  done.  For  that  case,  Gp  vs.  angle  was  included  in  the 
averaging.  We  show  the  new  values  in  Table  1  where  we  assume  5  dB  gain  for  surface 
antennas  and  0  dB  gain  for  the  FFB  antenna. 
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2500  KM 


Figure  8.  Probability  of  Correct  Message  Receipt  Vs.  Broadcast  Time  FFB-FFB;  PfGtG 


TBRD-5 


Figure  9.  Probability  of  Correct  Message  Receipt  Vs.  Power  Factor  P^G,.;  FFB-FFB;  Broadcast 


2500  KM 


Figure  10.  Probability  of  Correct  Message  Receipt  Vs.  Power  Factor  P+.G+Gr'»  FFB-FFB,  Broadcast  Time 


RVERRGED ,  TBRD=5 


Probability  of  Correct  Message  Receipt  Vs.  Power  Factor 


RVERRGED 


Probability  of  Correct  Message  Receipt  Vs.  Power  Factor  PtG+;  FFB-GRD;  Broadcast 


Probability  of  Correct  Message  Receipt  Vs.  Power  Factor  PtGtG  ;  GRD-GRD;  Broadcast  Time 


400  K 


Probability  of  Correct  Message  Receipt  Vs.  Power  Factor  P+GJ-i  ;  GRD-GRD;  Broadcast 


Table  1.  Radiated  Power  as  a  Function  of  Link  Connectivity 
and  Broadcast  Time  Using  Improved  Model 


Type 

Radiated  Power  (watts) 

PMR  = 

80% 

90% 

95% 

FFB-FFB 

5  min 

1100 

3200 

10  min 

475 

675 

910 

FFB-Surface 

5  min 

700 

1640 

3200 

10  min 

240 

420 

800 

Surface-Surface 

5  min 

220 

450 

725 

10  min 

70 

110 

172 

For  the  case  of  the  buoy,  we  are  interested  in  lower  powers  and  longer 
broadcast  times.  We  therefore  include  in  Figures  15  and  16  the  cases  of  1400  km 
GRD-GRD  for  30  and  60  minutes  respectively.  Using  5  dB  gain  for  the  receiving 
radiated  power  for  various  probabilities  of  message  receipt  for  various  broadcast 
times  in  Table  2. 

We  find  power  requirements  in  Tables  1  and  2  to  be  higher  than  that 

previously  posed  only  at  the  lower  values  of  power  as  in  the  10  min  FFB-FFB  case 

and  in  the  30  and  60  min  cases  for  the  buoy.  As  previously  mentioned,  at  values 

of  P.  G.G  >1000  w  the  old  model  based  upon  underdense  theory  is  adequate.  It  may 
t  t  r  “ 

even  be  conservative  since  at  large  values  of  T  and  lower  values  of  SQ,  the  fact 
that  overdense  bursts  have  longer  time  constants  can  overcome  their  relative  rarity 
by  their  higher  useful  life.  In  fact,  the  prior  time  T  used  in  this  analysis 
should  be  limited  to  a  few  seconds  because  trails  more  than  several  seconds  old  are 
unreliable  because  of  the  multipath  effects  due  to  wind-shear  distortion  of  the 
trail.  This  did  not  impact  our  analysis  since  the  largest  trail  time  we  use  is 
T+At  =  9At  =  1.097  s. 


TBRD=30 


Probability  of  Correct  Message  Receipt  Vs.  Power  Factor  P+G+G,,;  GRD-GRD;  Broadcast 


400  KM,  TBRQ-60 


Probability  of  Correct  Message  Receipt  Vs.  Power  Factor  Pf6fG  :  GRD-GRD;  Broadcast  Time 
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APPENDIX  A 

VHF(MB)  Link  Performance  Model 


The  VHF(MB)  model  Is  based  on  test  system  data  described  In  Refs.  1 
and  2  for  surface-to-surface  links.  The  effects  of  elevating  one  or  both 
antennas  is  estimated  from  geometrical  considerations.  These  aspects  of  the 
model  are  easily  modified  as  more  data  become  available. 

The  model  output  Is  the  probability  of  message  receipt  (PMR)  within 
a  continuous  broadcast  time  t.  The  message  may  be  In  a  single  packet  of  time- 
length  At  sec.  The  message  is  not  considered  received  unless  the  meteor  burst 
trail  is  of  sufficient  cross  section  and  duration  to  reflect  the  entire 
message.  For  a  message  consisting  of  N  packets,  the  probability  of  message 
receipt  is  defined  as  the  probability  of  receiving  all  N  packets  In  any  order 
within  the  broadcast  time. 

The  model  Is  based  on  the  physics  of  underdense  reflection  of  radio 
signals  as  described  in  Ref.  3.  An  underdense  trail's  diffusive  time  spread 
causes  the  received  signal  strength  to  decay  exponentially  with  characteristic 
time 

( sec) 
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where  x  «  wave  length  (m) 

2 

D  *  diffusion  coefficient  (m  /sec). 

In  the  model,  D  *  8  m  /sec  as  used  In  Ref.  1  and  +  ,  the  half  angle 
of  scatter  at  the  trail.  Is  calculated  for  a  meteor  burst  at  100  km  altitude 
midway  between  antennas.  Antenna  heights  are  zero  for  surface  and  21.3  km  for 
balloon.  For  surface-balloon  links,  the  calculation  uses  half  the  height  of 
balloon. 


A-l 


Although  useable  meteor  trails  occur  over  a  large  portion  of  the 

3 

sky,  the  most  likely  occurrences  are  near  the  midpoint  so  this  point  is  used 
for  scaling  purposes. 

When  at<<T,  a  meteor  trail  with  electron  line  density,  q,  high 
enough  to  scatter  the  required  signal  strength  at  the  receiver,  will  reflect 
the  entire  message.  We  designate  the  average  rate  of  occurrence  of  one  or 
more  bursts  with  this  value  q  or  higher  as  Mq  ( bursts/ml n) . 

3 

It  is  known  that  this  burst  rate  is  proportional  to  1/q.  For  under- 
dense  bursts,  the  scattered  signal  strength  is  proportional  to  q.  Reference  1 
outlined  a  method  of  scaling  from  their  test  system  to  any  other  test  system 
operating  at  optimum  ground  ranges.  For  a  required  receiver  threshold  Tr 
(watts)  which  sets  the  required  signal  strength  Tr>  and  given  values  of 
transmitter  gain,  Gt,  above  an  isotrope,  transmitter  radiated  power,  Pt,  and 
receiver  gain,  Gr,  above  an  Isotrope,  the  required  value  of  q  Is  proportional 
to 
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from  the  transmission  equation.  *  *  This  equation  is  shown  later  in  the 
text.  The  gains  correspond  to  directions  from  antenna  to  meteor  burst  which 
are  taken  at  the  midpoint  of  the  range.  Therefore  the  rate  of  useful  meteors. 


For  the  test  system,1 
x  *  6.383m 


and 


Wr 


A-2 


The  values  of  MQ  for  the  test  system  ranged  from  a  dally  minimum  of 
19/hour  to  25/hour  depending  on  season.  The  dally  median  value  ranged  from  50 
to  70/hour.  We  have  chosen  a  conservative  value  of  40/hour  =  2/3  /  minute  for 
the  model . 
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It  is  known  *  *  that  surface-to-surface  meteor  burst  performance  is 
best  at  midranges  but  degrades  at  larger  Tanges  because  of  earth  curvature  and 
at  shorter  ranges  because  of  geometric  and  depolarization  effects.  The  useful 
meteor  rate  calculated  is  therefore  modified  according  to  range  and  link 
geometry. 


A  linear  fall  off  of  useable  meteors  with  range  is  used  to  model  the 
decrease  In  the  volume  of  useable  meteor  trails  as  the  bulge  of  the  earth 
shadows  the  line  of  sight.  If  one  degree  elevation  is  used  as  the  limit  below 
which  signal  strength  drops  quickly,  the  groundrange  at  which  a  meteor  trail 
altitude  of  100  km  reaches  this  elevation  Is  1000  km  for  an  antenna  on  the 
earth's  surface.  We  therefore  choose  2000  km  range  between  surface  antennas 

as  a  cut  off.  The  actual  horizon  (0°  elevation)  is  at  1130  km.  According  to 

4  2 

Oetting,  a  fit  to  data  from  Boeing  Alaska  tests  shows  the  meteor  rate  at 

large  ranges  to  be  constant  for  770<R<1280  km  and  to  fall  off  for  R>1280  km. 

We  therefore  linearly  degrade  the  meteor  rate  to  zero  at  2000  km  from  its 

value  at  1280  km.  For  R<770  km,  we  use  Oetting' s  scaling  from  the  tests  from 

770  km  down  to  the  range  of  200  km  where  the  ground  wave  Is  of  sufficient 

strength  to  preclude  the  need  of  a  meteor  burst.  This  scaling  from  data  shows 

the  meteor  rate  reduced  a  factor  480/770  for  R<480  km  and  rising  as  R/770  to 

the  constant  value  for  480<R<770  km. 

For  the  case  of  a  surface  to  balloon  link,  the  ground  wave  can  be  of 
sufficient  strength  for  a  much  larger  range.  The  llne-of-slght  from  a  balloon 
at  21.3  km  altitude  grazes  the  earth  at  a  range  of  520  km.  Without  the  appro¬ 
priate  formulas  for  signal  strength  vs.  range  for  the  various  geometries,  we 
assume  that  200  km  of  further  ground  range  can  be  added  to  this  so  that  meteor 

bursts  are  not  needed  out  to  a  range  of  720  km.  Beyond  this  range,  the  effec- 
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tlve  meteor  burst  rate  will  be  taken  as  constant  as  in  the  surface-surface 
case  until  the  earth's  bulge  effect  sets  in  wherein  a  linear  fall  to  zero  at 
some  maximum  range  is  effected.  The  line-of-sight  from  a  balloon  to  a  meteor 
trail  grazes  the  earth  at  a  combined  range  of  520+1130  *  1650  km  from  balloon 
to  trail.  Since  we  choose  1000  km  as  the  maximum  range  from  a  surface  antenna 
to  the  meteor  trail,  the  total  cut  off  range  for  surface  to  balloon  is 
1650+1000  *  2650  km. 

The  determination  of  the  range  where  the  fall  off  in  useful  meteor 
trails  begins  requires  a  knowledge  of  the  distribution  of  the  trails  over  the 
sky  at  100  km  altitude.  The  only  case  where  this  distribution  is  given  in  the 
literature  is  for  surface-to-surface  at  1000  km  range.  This  shows  the  bulk  of 
the  useful  trails  to  be  near  the  midpoint  of  the  sky  between  the  antennas. 
This  should  be  also  true  at  greater  ranges.  Therefore  the  fall -off  with  range 
should  be  small  until  the  final  area  becomes  shadowed.  If  this  area  is  simi¬ 
lar  in  range  extent  to  that  of  the  surface-surface  case,  the  range  over  which 
the  fall -off  occurs  should  be  similar.  Me  therefore  begin  the  linear  fall -off 
720  km  in  from  the  maximum  range  of  2650  km  at  1930  km.  Similarity  in  range 
extent  is  consistent  with  using  the  same  constant  meteor  rate  for  both  surface- 
surface  and  surface-balloon  geometries. 

For  the  balloon-balloon  case  we  consider  ground  wave  operable  to  a 
range  of  520+520+200=1240  km.  The  cut  off  for  meteor  bursts  will  be  at  the 
maximum  range  1650+1650=3300  km.  We  take  the  meteor  burst  rate  as  the  same 
constant  out  to  a  range  3300-720=2580  km  where  the  linear  decline  begins. 
Figure  A-l  illustrates  the  modeled  range  dependence  of  MQ  for  the  three  cases. 

For  ranges  less  than  shown  In  the  figure  where  meteor  bursts  are  not 
needed,  the  model  returns  a  probability  of  message  receipt  of  unity.  For 
ranges  greater  than  shown  In  the  figure,  the  model  returns  a  probability  of 
zero.  The  geometries  for  the  minimum  and  maximum  ranges  are  shown  in  Figure 
A-2. 
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Figure  A-l.  Reference  Useful  Meteor  Burst  Rate 
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The  useable  meteor  rate  for  the  system  of  interest  is  obtained  by 


In  using  the  model  in  this  paper,  the  systems  considered  had  x  =>  6.383m  and  a 

receiver  detection  threshold*  of  T_  =  10'*^  /  data  rate  \  (  A  watts 

_3  r  \40oo  bits/ sec/  \OS3/ 

for  a  bit  error  rate  of  10_J  for  coherent  PSK. 

The  useable  meteor  rate  therefore  scales  as  M  =  Mq  Gr  for 

x  *  6.383  m  and  data  rate  *  4000  bits/ sec. 

The  scaling  of  M  and  the  exponential  time  dependence  of  the  signal 
which  will  be  used  in  deriving  the  PMR  output  of  this  model  is  valid  only  for 
underdense  meteor  burst  trails.  If  underdense,  the  signal  level  has  a  value 
and  time  dependence  which  varies  as 


-15 

where  r  =  2.8-10  m  =  classical  electron  radius, 
e 

3 

If  overdense  the  variation  is 


1  ft  In  /4reqT\ 
/4ir  t  y  t  J 


1/4 


4r  qx 

which  rises  to  a  maximum  at  t  *  — z —  of  value 


1/4 


and  drops  to  zero  at  t  =  4  r.qt  at  which  time  the  trail  is  completely  under- 

14  e 

dense.  If  4rgq<l  or  q<10  /m,  the  underdense  description  should  be  adequate 

and  conservative.  If  the  power  factor,  P*G*G_,  Is  so  low  that  a  large 
14  t  t  r 

q>10  /m  is  required,  the  model  will  have  to  be  extended  to  include  overdense 

reflections.  Further  analysis  is  required  to  determine  the  limit  of  applica¬ 
bility  accurately  in  terms  of  system  parameters.  An  estimate  can  be  made  from 

12  3 

the  transmission  equation  *  * 


,3_2_  2 

t  t  r  x  q  r 

- T  e 

16* 


_ sin2a _ 

RtRr^Rt+Rr^1“cos2e  sin2*V 


where  o,  the  angle  between  the  electric  vector  at  the  meteor  trail  and  Rr,  the 
direction  vector  from  meteor  trail  to  receiver,  is  a  measure  of  depolarization 
loss.  The  distance  from  transmitter  to  meteor  trail  is  R^.  and  s  is  the  angle 
between  the  trail  and  the  plane  of  Rt  and  Rr.  The  factor  in  brackets  dimin¬ 
ishes  only  slightly  with  range  since  a  and  *  get  larger  with  range.  If  B=0, 

g  3 

the  bracket  is  «/*10  km  at  the  limit  of  large  range  ^2000  km  for  which 


PtGtGr  = 


60  watts 

2  1 
q  r. 
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for  x  *  6.383m  and  Tp  =  10  watts.  The  transition  point  qrfi  *  1/4  occurs 

for  P+G..G..  =  1  kilowatt  at  2000  km  where  the  earth  bulge  has  no  usable  meteors 

^  ^  *  -8-3 

anyway.  At  1200  km  range,  the  bracket  is  ^2*10  km  and  the  transition 

point  occurs  for  P^G^Gp  *  500  watts.  Until  more  analysis  and  data  comparison 

is  done  we  should  consider  PtGtGr  =  500  watts  as  a  lower  limit  for  validity  of 

the  model  for  long  ranges.  At  shorter  ranges,  the  model  is  valid  at  lower 

power  factors. 

We  now  derive  the  formulas  for  the  PMR  from  the  meteor  burst  para¬ 
meters  M  and  t  and  the  system  parameters  At,  N  and  t.  In  the  derivation  we 
will  redefine  M  as  MQ  so  as  to  define  a  new  M  *  M0e"At/^T* 
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Consider  a  message  of  time  length  at.  For  any  set  of  parameters 
describing  the  communication  system  there  Is  a  required  minimum  electron  line 
density  of  the  meteor  trail  for  the  successful  receipt  of  a  short  message 
(at*0).  The  reflecting  ability  of  the  trails  decreases  exponentially  with 
time  with  a  time  constant  t.  If  a  message  Is  not  considered  received  unless 
its  entire  length  Is  received  continuously,  the  minimum  strength  of  the  meteor 
trail  must  be  increased  a  factor  eAt/^T  for  a  meteor  burst  occurring  just 
before  the  message.  If  the  message  Is  continuously  repeated  for  a  time  t  and 
the  message  is  considered  received  if  any  length  at  or  better  within  the 
message  train  Is  received  during  the  transmission,  the  probability  of  receipt 
is  calculated  as  follows. 

The  average  rate  of  occurrence  of  meteor  bursts  with  a  specified 
strength  or  better  Is  Inversely  proportional  to  strength.  Therefore  if  M  Is 
the  rate  of  required  meteor  bursts  for  at=0,  MaM{Je  '  is  the  rate  for  at^O. 

The  probability  that  a  useful  meteor  will  occur  in  the  short  time  dt 
is  Mdt.  The  probability  that  one  or  more  useful  meteors  will  occur  in  a 
longer  time  t  is  l-P(t)  where  P( t)  is  the  probability  that  no  useful  meteor 
occurred.  P{t)  is  calculated  as  the  product  of  the  probabilities  of  N 
independent  events.  Each  event  Is  the  nonoccurrence  of  a  useful  meteor  In  the 
time  interval  t/N  with  event  probability  (1  -  Mt/N). 


Suppose  Instead  that  the  message  Is  sent  only  once.  Reception  requires  a 
meteor  burst  of  required  strength  or  better  just  before  the  message  or  one 
stronger  by  edt/,T  a  time  dt  earlier  and  as  prior  time  goes  on  the  required 
strength  Increases  exponentially.  Suppose  prior  to  some  time  t  before  the 
message  began  there  were  no  meteors  of  the  required  strength.  We  can  then 
break  the  Interval  as  before  but  in  this  case 


Here  the  probability  of  nonoccurrence  approaches  unity  as  the  position  of  the 
time  interval  recedes  many  time  constants  prior  to  the  message. 


In  P 


N 

1  x 
1*1  x 

In  V 


N 

11m  z 

In 

N—  1*1 

fT1)  x 

Mt  vl\ 

-Mt 

r  J  l  \ 

TTX  )  3 

TT 

x  +  0  lr) 

.  *-xN+1 
~TT~ 

1 1m  -Mt 

e‘t/Nt  (l-e_t/,l 

A- 

x  *  e 


-t/Nx 


*  -Mt  (1-e  ^T) 


F(t)  =  e"MT(1’e  t/T) 


which  becomes  independent  of  t  for  large  t/x.  Therefore  the  probability  of 

M 

receiving  the  message  is  P  *  1-P(»)  *  1  -  e  T. 

If  the  message  Is  sent  J  times  In  a  string  with  repetition  period  T 
*  Mat  the  probability  of  not  receiving  the  first  and  second  transmissions  of 
the  message  Is  the  probability  of  not  receiving  the  first  times  the  probabil¬ 
ity  of  not  receiving  the  second  if  the  first  wasn't  received.  If  the  first 
wasn't  received,  no  burst  of  adequate  or  higher  strength  occurred.  Therefore, 
certainly  no  burst  of  strength  adequate  or  higher  for  the  second  occurred. 


Therefore,  no  meteors  occurred  of  the  required  strength  prior  to  T  before  the 
second  transmission.  The  probability  of  not  receiving  either  Is  then 


P(l,2)  *  P(-)  P{T)  =  e"Mt(2_e  T/t) 

If  we  figure  the  total  transmission  time  In  complete  periods,  t  *  JT 
■  JNat,  the  probability  of  not  receiving  any  of  the  J  transmissions  Is 

F(J)  =F(-)  F(T)(J-1)  *  e“MTe"T/T  e"MtJ  {1'e~T/t)  . 


At  1  2  J 


If  two  messages  separated  by  (p-l)at  spaces  are  continuously  sent  In  a 
packet  of  length  Nat  (l^p<N-l),  the  probability  of  not  receiving  both  messages 
if  J  packets  are  sent  Is 


P(A  or  B)  *  P(J)  +  P(J)  -  P ( A  and  B)  where  P(A  and  B),  the  probability  of 
receiving  neither  message  Is 

P(-)  P(pat)J  P[(N-p)at)]J‘1 


=  e-MTe"(N'pUtA  e-MTJ(l-e‘pAt/T)  e-MtJ  (l-e'(N'p)At/t) 


Since  a  broadcast  time,  t,  *  JNat,  we  see  that 


MTa-(N-p)at/t 

P  (J  and  TD  *  e  Mte  P_  Pu  _ 

P  M-p 


where  P  ,  e-M.t  fl-e-|,4t/T]  /Met 
P 


^e-IN-plPt/r 


At  this  point,  we  will  Ignore  the  multiplying  factor,  em  , 

since  Mt<<1  and  for  later  simplicity  in  nomenclature 

so  P(A  and  8)  3  pppN_p» 

P(A  or  B)  =  2Pn  -  pppN_p»  approximating  P(J)  3  PM. 

We  then  have  for  one  message  sent  continuously  with  no  spacing  (N*l)  the 
probability  of  receipt, 

PMR  3  1-Pr 

For  two  messages  sent  continuously  with  no  spacing  (Ns2),  the  probability 
of  receipt  of  both  Is 


PMR  3  1  -  P(A  or  B)  3  1  -  2P2  + 


If  three  messages  are  sent  In  a  packet  Nat  with  separations  of  (p-l)at 
and  (q-l)at 


A  B 


C  A  B 


A  B 


|  pAt  |  qAt  | 


the  probability  of  not  receiving  the  three. 


P(A  or  B  or  C)  *  P(A  or  6)  +  P(C)  -  P(CA  or  B]  and  C) 

Since  P(A  or  S)  »  P( A)  +  P(B)  -  P(A  and  B), 

PC CA  or  B]  and  C)  »  PC A  and  0  +  P(B  and  C)  -  P(A  and  B  and  C)  , 

PC  A  or  B  or  C)  ■  PC  A)  +  PC  B)  +  PC  C)  -  P(A  and  B)  -  PC A  and  0 

-  P(B  and  C)  +  PC  A  and  B  and  C) 


3PN  -  PpPN-p 


Pp+q  PN-p-q  "  PqPN-q  +  PpPqPN-p-q 


If  the  three  are  sent  continuously  with  no  spacing  (N»3),  the 
probability  of  receiving  all  three 


PMR  *  1-3P3  +  3P ]PZ  -  PL3 


Continuing  the  analysis  to  N  messages. 


PMR  *  1  -  NPn  +  z  PlAj  and  A2)  -  £  P(A1  and  A?  and  A^  + 

...  +  (-l)m  £  P(A1  and  A2  and  ...  and  AJ||)  ...  +  (-1)NP1N 

where  £  P(A^  and  Ag  and  ...  and  A^)  Is  the  sum  over  all  combinations  of  m 
messages  In  N  slots  of  the  probabilities  of  the  combination  not  being  received. 


For  N  3  4 

£  PCAj  and  Ag)  3  4pip3  for  the  combinations 


2 

plus  2P«  for  the  combinations 


Similarly  z  P(A.  and  1L  and  A-)  *  4P^P2  for  the  combinations 


liinimBliiilmBimimBiifiiilil 


so  P( 4)  *  1-4P4  +  4PxP3  +  2Pg  -  4?^  +  P^ 


The  next  four  formulas  of  PMR  are 


N=5,  PMR 


N=6 ,  PMR 


N*7,  PMR 


1  -  SP5  *  SPjP4  ♦  SP2P3  -  SP12P3  -  SfJ*  ♦  , 

5Pi3P2  - 

1  -  *s  +  6P1PS  *  6P2P4  +  3P32  -  6Pl2p4  -  12P1P2P3 
-  2P23  +  6P13P3  +  9P1ZPZ2  -  6P14P2  +  Pj6 


1  -  7P?  +  7PjPg  +  7P2P5  +  7P3P4  -  7P1£P5  -  14P1P2P4 

-  7PxP32  -  7P22P3  +  7Px3P4  +  21P12P2P3  +  7 P^3 

-  7Pl4p3  -  l4PlV  +  7p15p2  -  Pl7 

1  -  8P.  +  8P7P.  +  8PCP_  +  8P-P,  +  4P  2  -  8P,P  2 

8  71  62  53  4  61 

-  lepgPjjPj  -  16p4p3Pj  -  sp4p22  -  sp32p2  +  sp^3 

+  24P4P2P12  +  12P32P12  +  24P3P22P1  +  2P?4  -  8P4Pj4 

-  32?^?^  -  16P23P12  +  SPjPj5  +  20P22P14  -  BPgPj6  +  P^ 


Finally,  In  the  event  that  more  than  one  relay  Is  transmitting  an  N 
packet  message,  the  PMR  Increases  not  only  because  of  multiple  chances  to 
receive  the  whole  message  but  also  the  possibility  of  receiving  some  packets 
from  one  relay  and  the  rest  from  other  relays.  As  before,  the  PMR  Is  given  In 
terms  of  sums  of  probabilities 


N*8,  PMR 


P  (A1  and  Ag  and  ...  and  A^)  *  [J  Pj 

J 

where  the  values  of  j  In  the  product  depend  on  the  positions  of  the  packets  In 
the  N  packet  message.  For  L  relays  transmitting  this  becomes 


L  L 

n  n  p,  (i)  - n  n  p,m. 

1*1  4  J  j  1-1  J 


The  PMR  Is  thus  calculated  as  before  for  one  link  with  the  sole  change  being 
Pj  replaced  by 


77  P<(1)  where  P.i(1)  Is  calculated  for  the  ith  link. 
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APPENDIX  B 
IMPLEMENTATION 


The  implementation  of  the  overdense  formulas  into  the  VHF(hffi)  Link 
Performance  Model  requires  the  redefinition  of  repeated  quantities  for  minimum 
running  time.  For  the  same  reason  the  flow  diagram  will  appear  repetitious.  For 
example,  the  overdense  formulas  are  shown  separately  for  the  cases  q-pl  and  q-j>l. 

For  the  q-]=l  case  we  require 

P-  (T-t  ,  At+t  )  in  which  case  our  defined 

O  O  0 


from  the  definition  of  t„  and  the  fact  that 

o 


At  is  At+t  for  this  case.  In  the  definition  of  A  =  — —  ,  the  sum,  T+At,  is 

0  D-B  qTT 

unchanged.  We  define  D  i  At/x  so  that  qQ  *  e  so  the  first  check  on  underdense 

or  overdense  is  D>B  in  the  flow  chart  on  Figure  B-l . 


Because  the  formulas  will  be  used  for  the  special  cases  T  *  PAt,  1 <P<N ; 
we  define  TT  =  (P+1)D  so  A  =  TT/qT.  In  the  underdense  case  D<B,  the  check  to  see 
if  overdense  formulas  are  not  also  required,  t^T.becomes  B>JT  by  substitution. 

If  the  check  is  positive,  the  underdense  formula  of  reference  1  is  used.  This  is 


p-  .  /V  NAt  Qp 
u 

where  Qp  =  (l-e~T//x)/q0.  The  term  is  the  multiplier  of  the  probability  of  no 
useful  meteor  trail  for  repetitions  where  t  is  the  broadcast  time.  Since  the 
quantity  that  differentiates  the  various  situations  is  Qp,  this  is  calculated  first. 


If  the  check  is  negative,  the  fact  that  the  probability  is  the  product 


n —  .  -V  NAt 


P»** 


and  the  overdense  PQ  is  accounted  for  by  adding  l/qQ-l  to  the  Qp  calculated  for  Pc 


In  this  case  the  check  for  the  situation  of  Figure  4  or  5  which  is  —  1<12 
becomes  B  >q2«  If  negative,  the  situation  of  Figure  4  applies.  The  next  check 

1^2  written  TT  <q2.  If  positive  we  must  obtain  qy  from  the  transcendental 
equation 

4  4(T+At)/ (xqT) 
qT  =  S0  e 

which  in  this  case  is  qy  =  e4(A_B)  =  e  (^T  )  ,  The  first  time  this  is 

encountered,  the  solution  should  be  near  unity  so  qy  is  initialized  at  unity. 

For  subsequent  entries  at  higher  values  of  P,  the  previously  obtained  qy  initializes 
the  next  search  for  the  solution. 

For  TT>q2,  we  have  entered  the  situation  of  Figure  5  and  must  obtain 
the  quantity  qM  from  the  transcendental  equation 

n  -  n  a4(TT-qM) 
qM  -  q2e  SM'  . 

Since  qM  =  q2  at  the  check  point  TT=q2  we  initialize  qM  at  q2- 

For  B_>q2»  we  are  in  the  situation  of  Figure  5  from  the  beginning.  The 
integration  is  between  q3=l  and  the  qM  found  as  before. 

For  q^l,  we  have  the  overdense  case.  If  D<q2,  we  must  redefine 
BB  =  (At/q^t)  =  (D/q-j )  where  q^  must  be  obtained  from  the  transcendental  equation 

a  *  e4(BB-B) 
q!  e 

We  use  qQ  as  the  initial  guess.  We  then  do  the  check  TT^.  If  positive,  we  find 
A  =  TT/qT  from  the  same  equation  as  before  but  this  time  the  initial  guess  can  not 
be  qT=l  since  we  know  q-|iqyiq2*  We  start  the  search  at  qy®q-|  and  start  subsequent 
searches  at  the  previously  obtained  qy.  For  TT>q2  we  obtain  q^  as  before. 

For  D>q2  we  must  find  the  starting  value  q3  from  the  equation 
4(D-q3) 

q3  =  q2e 


where  we  use  q3sq2  for  the  initial  guess. 


FIGURE  B1 
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Figure  B2.  Program  Flow  Chart,  Overdense. 
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